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Abstract 

Some direct relations between soliton solutions of integrable hierarchies and thermody- 
namical quantities of the Coulomb plasmas on the plane are revealed. We find that certain 
soliton solutions of the Kadomtsev-Petviashvili (KP) and B-type KP (BKP) hierarchies 
describe two-dimensional one or two component plasmas at special boundary conditions 
and fixed temperatures. It is shown that different reductions of integrable hierarchies 
describe one (two) component plasmas or dipole gases on one-dimensional submanifolds 
embedded in the two-dimensional space. We demonstrate application of the methods of 
soliton theory to statistical mechanics of such systems. 

Keywords: Coulomb plasmas; Integrable hierarchies; Tau-functions; Solitons; Fermion sys- 
tems; Dipole gases 

1. Introduction 

Recently we have shown that the grand partition functions of some one-dimensional lattice 
gas models or equivalent to them partition functions of some Ising chains, coincide with the 
A^-soliton tau-functions of various hierarchies of integrable nonlinear evolution equations. The 
present paper is the third one in the series and it comprises a detailed comparison of exactly 
solvable Coulomb plasma models on lattices with integrable equations. In extension of our pre- 
vious considerations we discuss statistical mechanics of the Coulomb (logarithmic interaction) 
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gases on intrinsic two-dimensional geometric figures and various one- dimensional submanifolds 
of the plane. In this way we do not merely reinterpret previously known results but also 

reveal a number of new exactly solvable models. It is natural to expect that the connection 
with integrable equations gives a clue to the classification of solvable (at fixed temperatures) 
plasma models. 

A classical Coulomb plasma is a system of charged particles interacting through the Coulomb 
potential. In the Euclidean space R"' the Coulomb potential is defined as a solution to the 
Poisson equation 

AV{r,r') = -u;J{r-r'), u;^ = ^f^, r e (1) 

T{n/2) 

with V{r, r') obeying certain boundary conditions. For the two-dimensional plane without 
boundaries, equation (|I]) has the following solution 

V{z,z') = - \n\z- z'\, (2) 

where z = x + ly and ct;2 = 27r. A system of particles forms a stable plasma if the opposite 
valued charges do not recombine with each other forming a gas of neutral molecules. 

In general the location of plasma is constrained to some domain in which case there is a 
nontrivial interaction with the boundaries of this domain. In particular, the normal component 
of the electric field 

S = -VV (3) 
should vanish on the surface of an ideal dielectric 

£n = 0, (4) 

while the tangent component of this field vanishes on the surface of an ideal conductor 

£t = 0. (5) 

A useful way of solving the Poisson equation with such boundary conditions is provided by 
the method of images. In the present work we consider systems where every charge has either 
a finite number of images created by the boundaries or boundary conditions create periodic 
lattices of images. In the case of a finite number of images the solution to (|l]) is a finite sum of 
logarithmic potentials (||) created by a charge and its images. 

The energy of such systems of particles has the following form 

En= E Q^Q,v{^^,z,)+ E QM^i)+ E QM^^) (6) 

l<j<j<Ar l<i<N l<i<N 



where Zi = Xi + lyi and Qi denote the coordinate and the charge of the i-th particle on the 
plane. 

The first term in (^) is the energy of interaction between different charges. The second term 
is the sum of self-energies. This term is constant at certain boundary conditions, e.g. when the 
plasma is homogeneous or if it is restricted to a line with transverse boundary conditions. In 
general it arises as the energy of interaction between the charge and its own images. The third 
term describes an interaction of charges with external fields. 

The grand partition function of a system of particles of s different species is 

^ — 2^ ■ ■ ■ 2^ j r^ni...ns 

= y ... y / e-r-f^'H-.+/^i"i+-+M."«(ifi, (7) 

where fii, . . . , fig denote chemical potentials (Cs = e^" are the fugacities) and Q is the integration 
measure over the configuration space occupied by particles. In (|^ we have introduced the 
dimensionless inverse temperature 

r = /3g^ p = i/kT (8) 

and the dimensionless Hamiltonian 

Hn^...n, = \H mjy{.Zi, Zj) + y q}v{zi) + y qi<p{zi), (9) 

ijkj i i 

where qi = Qi/Q and (j){z) = ^{z)/Q. For one and two component plasmas it is convenient 
to choose Q = \Qi\, so that the dimensionless charges = 1 in the one component case and 
qi = ±1 in the two component case. 

The two-dimensional one and two component plasma models have been solved for a vari- 
ety of boundary conditions at the special value of the inverse temperature F = 2 (see, e.g. 

and references therein). In the case of two component plasma an exact solution exists 
due to the well known correspondence between this system and the free fermion point of the 
Thirring/Sine-Gordon model. In technical terms the solution is possible due to different de- 
terminant representations (the Cauchy determinant for the two component plasma and the 
Vandermonde determinant for the one component case). Such determinant representations al- 
low one to solve models of log-gases on a line with the transverse boundary conditions (e.g., 
see [|]). Note that the main part of the cited works deal with the continuous space models. 
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for an account of the lattice models, see, e.g. ^ 0]. The literature on the Coulomb gases 
is enormous, many statistical mechanics models have been related to them |Q, there is a re- 
lation to conformal field theories, etc. Still, the identification of Coulomb plasmas on lattices 
with some boundaries and of the famous multi-soliton systems has been missed in the previous 
investigations. 



2. Basic observations 

Let us consider a lattice version of (|^. We suppose that each type of particles can occupy 
only a discrete set of points in the complex plane. E.g., in the two component case the positive 
and negative charge particles occupy sublattices {z^} and {z-}. We denote the union of all 
sublattices as {z}. No more than one particle is allowed at each site. In this case, the integrals 
in ([^ are replaced by discrete sums over the lattice points and the whole partition function 



can be rewritten in the following form ||T0| 



G = T.^wilT.W{z, z') a{z)a{z') + ^ nj{z)a{z) , (10) 

where 

W {z, z') = -Vq{z)q{z')V{z, z'), w{z) = fi{z) - T {q\z)v{z) + g(z)0(z)) (11) 

and a{z) = or 1 is an occupation number of the site with the coordinate z. The variables 
q{z),^{z) are some functions of the lattice coordinates. For example, q{z±) = ±1, /i(-2±) = fj,± 
for the two component plasma. 

Now, let us write out the r-function of A^-soliton solutions of some integrable hierarchy in 



the Hirota form 11 



T^= E exp i E A,M^)a{z') + E 0{z)a{z) , (12) 

(7=0,1 \ Z^Z' {z} J 

where the variable z takes discrete values describing spectral characteristics of solitons. The 
function 6{z) parameterizes phases of solitons with the index z and Azz' is the phase shift 
acquired as a result of the scattering of solitons with the indices z and z' on each other. 

Evidently, the expressions (0) and (|T2|) have the same form. One just needs to make 
proper identifications between the phase shifts A^z' and the interaction potentials W{z, z'), and 
between the phases 9{z) and the function w{z). 



Concluding this section we present brief explanations on the structure of the expression 
([T2I) . One can track appearance of the r-function on the simplest example of the Korteweg-de 
Vries (KdV) equation which emerges as a compatibility condition of two linear problems: 

LiIj{x) = -tpxxix) + u{x)tlj{x) = Xip{x), (13) 

ipt{x,t) = B^{x,t), B = -Adl + 6u{x,t)d^, + 3u^{x,t). (14) 

Eigenvalues of the Schrodinger operator ([T3|) do not depend on "time" t: dX/dt = 0. As a result, 
the compatibility condition of ([13D and ([I^) takes the following operator form: Lf = [B,L], 
and yields the KdV equation 

Ut + Uxxx - QuUx = 0, 

describing the motion of shallow water. Stable traveling wave solutions of this equation are 
called solitons. The inverse scattering method allowed one to find an explicit expression for the 
general iV-soliton solution 

u{x,t) = -2dl\nTN{x,t), (15) 

where is the determinant of some N x N matrix called the tau-function. It can be written 
in the form 

Ati=0,l \l<i<j<N l<i<N J 

Oi = hx - + ef\ z, J = 1,2,..., AT. (17) 

The variable ki is related to the i-th eigenvalue of the operator L in the standard quantum 
mechanical meaning, Aj = —kf/4, and it parameterizes the amplitude of the i-th soliton which 
is proportional to kf. The scattering process is described by the time evolution from t —00 
to t — s> +00. In this picture /ki are the zero time phases of solitons and kf are equal to their 
velocities. 

The phase shifts Aij describe relative retardation of i-th soliton after the interaction with 
the j-th one. It is determined by the formula 

{ki ~\~ ^^j)'^ 

In the standard approach ki are real. However, formally one can take them as complex 
variables as well. Then it is seen from (0) that the plane coordinates Zi of some particular 
configuration of the Coulomb particles (see below) play the role of ki and correspond thus to 
complex amplitudes and velocities of solitons. This leads to complexification of the function 
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u{x, t) which requires somewhat different physical interpretation of the KdV equation. However, 
for complex ki there is no complete coincidence of (0) and (12) — the signs of modules in 
the Coulomb potential make G (|TUp real and there are no such signs in (|TB]) and, so, tn is 
complex. On the one hand, one may place the plasma onto the x or y axis in order to reach the 
coincidence and this one-dimensional situation was considered in On the other hand, one 
may try to identify KdV solitons with complex ki with a special system of interacting electric 
and magnetic charges ||^. We shall not consider the latter situation here, but rather concentrate 
upon the purely electric systems on the plane with real energy of interaction. 

There are generalizations of the expressions ([T7| ) -([T8| ) such that the corresponding u{x, t, . . .) 
satisfy higher order members of the KdV-hierarchy, sin-Gordon, Kadomtsev-Petviashvili (KP), 
Toda, and some other integrable equations |Tl|. Note that all these hierarchies can be derived 
within the free fermion formalism which is outlined in the concluding section of the present 
article. In the next section we consider in detail an identification of the KP hierarchy solitons 
and plasma particles. 



3. KP hierarchy 

The KP hierarchy which starts from the KP equation 

d'^u d f du g ^"^ _|_ d^u\ 
dy"^ dx \ dt dx dx^ J 

is obtained through a generalization of the KdV formalism. The KdV equation is an equation 

which describes the eigenvalue preserving deformations of a second order differential operator, 

while the KP equation describes that of a first order pseudo-differential operator: 

L{t,d)^{t,k) = kip{t,k), 

k) = Br.it, d)^it, k), 

dtn 

where t = {ti = x,t2 = y,t^ = t, . . .), d = and 

L{t, d) = d + ui{t)d-^ + u^{t)d~^ + . . . , 

B,,{t,d) = d- + Y.h^,{t)dK (19) 

j=0 

The KP hierarchy is a system of nonlinear differential equations for Uj{t) and bnj{t), resulting 
from the Zakharov-Shabat type compatibility conditions of linear equations given above 

dL _ dBrr. dB^ _ FR R 1 

— - [Bn, L\ , — -— - [Bn, Bm\ ■ 
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As in the KdV case, the potential u = Ui is expressed in terms of the tau function (|T5D and 
iV-sohton solution is presentable in the Hirota form (|l^) with the following parameterization 
of the soliton phases and phase shifts [|T^]: 

, , {az-a.z'){hz-h^,) 

oo 

e{z) = e^'\z) + Y.i<-i-bzr)tp, (20) 
p=i 

where tp is the p-th KP "time" and a^, are some arbitrary functions of z. 
If we take 

ttz = z = X + ly, hz = —z* = —X + ly, y >0 



then 



where 



Azz' = Wiz,z') = -2Viz,z'), (21) 



V{z, z') = — In 1^ — z'\ + In \z* — z'\ 

is the potential at the point z created by a positive unit charge particle placed at the point 
z' over the conducting surface occupying the y < region (since V{z, z') solves the Poisson 
equation with the tangent boundary condition (||) at ?/ = 0). Using the method of images one 
may say that this is an effective potential created by a positive charge at the point z' , '^z' > 0, 
and its image of opposite charge located at {z')*. Comparing (|2l|) with the original definition 
( prp one finds that the temperature (|]) is fixed and equals to 

r = 2. (22) 

Thus the reduction of KP A^-soliton solution described above corresponds to the Coulomb 
plasma in the upper half plane '^z > with metallic boundary along the x-axis. This situation 
is depicted in the Fig. 1. 

Note that the choice = z,bz = z* gives a similar situation — a plasma in the right half 
plane with the metallic boundary along the y axis (evidently this configuration is reached by 
the rotation of the coordinates z iz). 

Let us shift z ^ z + la, a real and take the limit a — >■ oo, i.e. take the plasma far away 
from the boundary. This leads to some divergences in the energy which can be removed by 
addition of an appropriate diverging constant to the initial Hamiltonian. As a result one gets 
the pure plasma system at the inverse temperature F = 2. This temperature corresponds to the 
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Figure 1: KP equation: one component two-dimensional plasma above an ideal conductor. 
Positive charges are shown as white squares while their negative images are shown as black 
squares. Interactions between different charges are shown by dashed lines, while the interactions 
between charges and their own images are shown by the solid lines. 

standard normalization in the random matrix theory [^]. Note that the normalization of the 
temperature in our previous papers was chosen as F = 1 since there we were discussing Ising 
chains without detailed comparison with the Coulomb systems which is a goal of the present 
work. 

The identification w{z) = 6{z) allows one to write the following expression for the zero-time 
phases O'-'^^z) in (|20|) 

e^^\z) = fi- r {In \z* - z\ + (f){z)), r = 2, (23) 

where the second term corresponds to the "charge-image" self-interaction energy, and the last 
term describes the potential of the field created by the neutralizing background of the density 
piz): 

0. (24) 



A0(z) = -27rp(z), 1^ 



y=0 



Contributions from the KP "times" 



Y,{zP-{zy)tp = -Tct>,^t{z) (25) 

p=0 



correspond to an external electric field. One can see that ( pSj) is a sum of the harmonic 
polynomials, p-th polynomial being a fundamental antisymmetric polynomial of the group 
of reflections (the Coxeter group) of 2p-gon. Since the Laplacian of this part is zero the 
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corresponding density of charges is zero, i.e. this field is generated by external distant charges. 
For instance, the contribution of the first time {z — z*)ti corresponds to the homogeneous 
electric field perpendicular to the boundary. For reality of the potential the time ti has to be 
purely imaginary. 

The system of distant charges Qi located at the points Wi above the conductor create the 
following electrostatic potential at the point z 



Since the external charges are far from the origin \z\ <S \wi\, we can expand the potential 
4>ext{z) in the Taylor series and get as KP times 



t ^-V ( ^ ^ 



p 



In this picture the KP times take imaginary values automatically. 

One may conclude that the general imaginary times evolution of a special system of KP 
hierarchy sohtons describes electrostatics of a plasma in a varying external electric field. 

It is worth of mentioning that sometimes the continuous limit is not physical below some 
temperature. Indeed, let us consider the partition function of particle gas at the inverse 
temperature F in a square domain {0 < x — ^z < L) x {0 < y — 'isz < L) 

Zn^I dxi [ dyi... [ dxN I dt/iv exp ( F V In -p^^ — ^ - F Vln l^j - | . 

Jo Jo Jo Jo \ i<3 i J 

Introducing new dimensionless variables Zi/L one gets the following expression for the thermo- 
dynamic potential or the pressure P in appropriate units 

„ IuZn (2-F)lnL 

P = — — h terms independent on L. 

We see that at the inverse temperature F = 2 the average pressure changes the sign. The 
system is inhomogcncous and the local pressure is different from the average one. It is clear 
that the local pressure changes the sign at the boundary and at the inverse temperatures F > 2 
the particles stick to the surface of the conductor. 

As a result, the continuous space version of the model in a domain with the conducting 
boundary does not make physical sense if the temperature is low enough. However, the situation 
is different if plasma is confined to a domain which does not touch the boundary (as in the 
situation depicted in the Fig. 1). E.g., one can avoid such an unphysical behavior in the 



discrete (lattice) version of the model. In the latter case the lattice spacing plays the role of 
the radius of a hard core repulsive interaction which prevents collapse of the system. Note that 
the identification with the KP soliton systems takes place exactly at the critical temperature 
r = 2. 

Closing this section let us note that the equation (|l]) and the boundary conditions (H), (||) 
hold for an appropriate conformal change of the variable z —>■ f{z). For instance, choosing 
soliton parameters in ( PPD as follows = z^, = — (z*)^, we obtain 



W{z,z') = 21n 



z^ - (z'Y 



{z*y - (z')^ 

corresponding to the interaction of two charged particles in the rectangular corner with metallic 
walls along the x and y axes. Higher degree monomial maps z ^ z"" put the plasma into the 
corner with the ir/n angle between the conducting walls. 
The exponential map 



exp(7rz/L), bz = — exp {ttz* / L) 



(26) 



generates the VT-potential 



W{z,z') = 21n 



sinh ^{z — z') 



sinh ^{z* — z') 

which describes the plasma in the strip ^z = (0, L) between two parallel conductors. 



(27) 



The choice 



results in 



ttz = exp{—TTx/L) 
W{x, x) = In 



bz = — exp(— 7r(a; + a)/L), 



sinh 



2 ■k{x-x') 
2L 



sinh ^ {x 



(28) 



X' — a) smn -^[x — x' -\- a) 
Solution of the Poisson equation with periodic boundary condition along the ?/-axis with the 
period 2L is given by the potential M 



V{z,z') 



In 



sinh^(z 



Therefore one can interpret (|2^) as the interaction energy of two neutral dipoles in the periodic 
background with the distance between charges in the molecule equal to a (the internal energy 
of dipoles is neglected). These dipoles all lie on the x-axis and have identical orientation. Since 
W = —TV, we see that the effective inverse temperature is twice smaller than in the previous 
cases, i.e. P = 1. 
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Let us choose parameters in (pO|) as specific double periodic functions 

2 I. 2 * 

Qz = sn z, Oz = — sn z , 



where sn z is the Jacobian elhptic function with the periods Lx,iLy. Then one gets the W- 
potential 



W{z,z') = 21n 



2 2 / 

sn z — sn 2; 



sn^ z* — sn^ 



2 In 



hiz - z')e,{z + z') 



ei{z* - z')ei{z* + z') 



(29) 



where 9i{z) is the Jacobi theta-function 



^i(z) = 2^(-l)"exp -7r(n + l/2) 



n=0 



Ly\ , 7r{2n + 1)2; 



sm 



Lx 



It vanishes when z lies on the boundary of a L^; x Ly rectangle, i.e. we have the plasma in a 
box with the conducting walls (F = 2). For small z, z' (or for the large size box L^,, Ly 00) 
one recovers plasma in the rectangular corner. 



4. KP hierarchy: the two component plasma case 

Let us consider the situation when lattice points {z} consist of two subsets {^i}, {-2} = {-2-} U 
{2+}. Choosing the following identification of parameters in ( pOD 

az^ = z^, bz_ = -z*_, = z^, b^^ = -z+, (30) 

we get a model of the two component plasma (the plasma consisting of two species of opposite 
charges) above an ideal conductor at the inverse temperature F = 2: 

W{z±,z'^) = 21n|2± -41 -21n|4 -41, 

W{z±,z'^) = 2 In 14 -2^1 -21n|2± - 4|, 

G{z±) = /i± - 21n \z± - z^l =F 2(/)(2:±). 

After the conformal transformations z — ^ z^,e^ , etc one gets the two-component plasma in the 
metallic rectangular corner, a strip, etc. 

Note that the bulk properties of the two component plasma are different from the ones of 
the one component systems. Actually, statistical mechanics of the one component plasma is not 
stable without metallic boundary or neutralizing background, since its particles tend toward 
boundaries repelling each other with long range (logarithmic) forces. In the two component case 
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screening is possible if the system is neutral and the temperature is high enough. Simple scaling 
analysis shows that the pure two component plasma undergoes a transition in the bulk at F = 2 
0. At this temperature an association of opposite charges into molecules (neutralization of 
plasma) takes place. Again, there is no such a transition in the lattice case because of the 
hard-core repulsion at the lattice spacing distance. 

It is well known that the two component homogeneous plasma without boundaries is a 
specific representation of the quantum Sine-Gordon or Thirring model \T3\. The inverse tem- 
perature (^) corresponds to the free fermion point, but the model is (in principle) integrable 
for any F. The Sine-Gordon model requires renormalization if the coupling constant exceeds a 
critical value. The lattice spacing (in the case of the lattice plasma) is equivalent to the intro- 
duction of a cutoff in the corresponding field theory. Thus it is natural to associate the case 
of two component lattice plasma above the metallic boundary with some discretized boundary 
Sine-Gordon or Thirring model. 



5. KP hierarchy: some reductions 

In this section we consider a number of one- dimensional reductions of the KP hierarchy and 
some self-similar soliton solutions. We describe here only the most popular integrable systems 
and do not cover all possible cases and their Coulomb gas interpretations. 

1) We begin with the reduction considered earlier in the literature [Q. In this case, plasma 
is restricted to the line y = Y. The IV-potential is translationally invariant and equals to 

+ [x — X Y 

Self-interaction energies of charges with their own images are constant and may be neglected. 

2) Reduction to the KdV hierarchy. In this case particles move along the vertical line a; = 
and the VT-potential is 

W^(y,y')=lnj^^. (31) 

[y + y) 

Now a non-trivial self- interaction term oc In |2?/| enters the definition of soliton phases. 

3) Discrete KdV hierarchy. The phase shifts have the form 

Taking z to be purely imaginary z = la we get the following result 



W{a,a') = In 



/ sin |(a — a') \ ^ 
\sin ^{a + a') 

12 



which corresponds to the plasma restricted to an arc with the center of the corresponding circle 
at the conductor's surface. The self-interaction energy is — ln|2sina|. It is sufficient to put 
z = e'° in ( piD in order to get this system from the KP soliton solutions. 

4) The reductions admitting both left and right moving solitons correspond to the plasma 
restricted to domains of disjoint parts. E.g., the Toda lattice case |jlT[ 

, , / e(z)z — e(z')z' \ ^ , , 
y 1 — e{z)e{z')zz' J 

corresponds for z = e^" to the two component plasma on the half circle with the center lying 
upon the conductor surface. The positive charge particles occupy the subsector a G [0,7r/2] 
and the negative charge particles are located in the subsector a G [7r/2,7r]. 

5) The Boussinesq equation corresponds to a one component plasma on disjoint halves 



of the hyperbola 3x = y + 1 situated above the conductor occupying the y < half plane: 



It would be interesting to find the equation whose soliton phase shifts are given by the 
elliptic functions (pOl) . 

Consider some specific lattice plasma configurations associated with self-similar soliton so- 
lutions of integrable equations. 

Let us take the KP hierarchy with parameters (|26|) and restrict the corresponding plasma 
to the line parallel to the conductor surfaces 

^z = Y. 

Applying this constraint, we get from (27) 



W{x - x') = - In sm^ — coth^ — - + cos^ — . 

The simplest self-similar reduction corresponds to the homogeneous one-dimensional lattice 
parallel to the x-axis. Then the n-th charge has the coordinates Zn = X + hn + lY , where X is 
some fixed constant and h is the lattice spacing. In this case soliton momenta and hz form 
one geometric progression (|26|) . The simplest KdV self-similar reduction corresponds to the 
case when the line is located at equal distances between parallel conductors Y = L/2: 

W{x - x') = In tanh^ ^i^-^) _ ^^^^ 
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M-periodic reduction corresponds to the situation when sohton momenta are composed 
from M distinct geometric progressions. In this case the lattice {z} consists of M homogeneous 
sublattices. In the plasma language it corresponds to the model where plasma moves on M 
distinct parallel lines between two conductors with the coordinates {z} = {lYp + Xp + nh,p = 
1, . . . , M, n = 0, 1, . . .}. If one sets Yp = L/2 then all these sublattices are situated upon the 



middle line and this case corresponds to the general self-similar KdV soliton potentials of |T6 

In a similar way one can place plasma upon M parallel lines with fixed x-coordinates, 
X = Xp, which are situated above the y < conductor. In this case self-similar lattices are 
described by restriction of y to the union of M geometric progressions {z} = {Xp+iYpg", q < 1}. 
However, such configurations are not safe from the collapse of particles on the walls. The KP 
self-similar systems are richer than the ones of the KdV or BKP (to be considered below) cases 
due to the presence of non-trivial translational parts Xp (or Yp) in the parameterization of the 
corresponding soliton spectral data az,bz- 



6. Correspondence with Ising lattices 

It is well known that lattice gas models are related to the Ising models [[l^. Indeed, substituting 

a{z) = {s{z) + l)/2, s{z) = ±l, 

into (|l3); we get a formula for the partition function of an Ising model up to some constant 
multiplicative factor 

g = T.^w{Ij:j (^^ ^i^H^') + E ^(^)^(^) I , (33) 

H V ^i^z' {z} J 



where 



are the exchange constants and 



J{z,z') = ^W{z,z') 



H{z) = ^w{z) + \ Y: W{z,z') (34) 

z' ,z'=^z 

denotes the external magnetic field. Here we have absorbed the temperature variable F into 
the definition of J{z,z') and H{z). 

Note that the Ising models emerging in this way look natural only in the reduced one- 
dimensional case, when there are no constraints upon the values of spins at different points. 
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In the two-dimensional picture with boundaries one has to assume that the configuration of 
spins satisfies some geometric constraints which do not have natural meaning similar to the 
one existing in the electrostatics. This leads to the exchange which depends not only on the 
distance between the spins but on the distance to the boundaries as well. 

Technically, it appears that the Ising representation of the lattice plasma grand partition 
function is more convenient for writing it in the determinant form. 



7. BKP hierarchy 

Consider the BKP equation: 



— I 9— - St^t^ + — - 30—— + 30—— + 60 I — I I - 5^ = 0. 



dtA dh dhdty dt\ dhdti dhdtl^ \dtJ I dtl 



This equation involves three independent variables ^1,^3,^5. It is the first equation of the B-type 
KP (BKP) hierarchy which is a reduction of the general KP hierarchy [|T^, [T^. This reduction 
is achieved by setting bno = in (|19[), which assumes, in turn, that even "times" evolution is 



absent, t2p = 0. The BKP tau-function is defined via the relation u{ti,t^, ■ ■ ■) = dlnrsKp/dti. 

A^-soliton solution of the BKP hierarchy has the same form (0) with the following phase 
shifts and soliton phases 

^ ^ _ (Qz - a'z'){bz - b,>){a^ - h^i){h^ - a^>) 

00 

e{z) = e^'\z) + E(ar^ + fer^)t2p-i, (35) 

p=l 

where t2p-i are the BKP "times" . 

If we take one component plasma and fix 

ttz = z = X + ly, bz = z* = X — ly, x > 0, y > 0, 

then Azz' = W{z,z') = -2V{z,z'), where 

V{z, z') = - In |z - ^'1 - In 1^* - ^'1 + In |z + /| + In |z* + z'\ 

is the total interaction energy between unit charge particles at the points z, Qz > 0, > 0, 
and z', '^z' > 0, 3?z' > 0, in a domain of the upper right quarter of the plane with an ideal 
dielectric boundary (cf. (H)) along the x-axis and an ideal conductor wall along the y-axis (cf. 
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Figure 2: BKP equation: a one component two-dimensional plasma in the corner between an 
ideal dielectric (the horizontal axis) and an ideal conductor (the vertical axis) . Positive charges 
are shown as white squares while their negative charge images are shown as the black squares. 

Interactions between different charges are shown by dashed lines, while the self-interactions 
between charges and their own images are shown by the solid lines. 
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(1^)). This situation is depicted in the Fig. 2. According to the inverse temperature is fixed 
and equals to F = 2. 

If the plasma is taken far away from the corner by appropriate translations, then one gets 
the pure plasma system at the effective temperature F = 2. Sliding along the y-axis to infinity 
one comes to the previously considered plasma associated with the KP equation. Sliding along 
the X-axis requires a renormalization of the zero energy level after which one gets a plasma 
above the surface of a dielectric. 

If we place charges upon the y = axis, then the dielectric boundary condition disappears 
and we have 

W{x,x') = In 



X — X 



X + x' ^ 

which corresponds to the plasma model induced by the KdV equation at the inverse temperature 
F = 4, which is twice higher than in the appropriate KP reduction case. 
Identification of the initial phase 9^^\z) in ( pOD is as follows 



e^°\z) = 21n| z* — z\ — 2 In \z* + 2;| — 2 In \2z\ + /i — 2(f){z), 
where the meaning of all terms is similar to that of (|23|), (p^). Contribution of the BKP "times" 

oo 

Y^iz'^-' + {z*yP-')t2p-, (36) 
p=l 

corresponds to electric field created by external charges. The potential (^) satisfies all the 
boundary conditions @), The p-th polynomial in the sum (|36|) is the fundamental anti- 
symmetric polynomial of the refiection group of the 2(2p — l)-gon. Again, the evolution of 
solitons under the BKP fiow corresponds to the evolution of plasma under a motion of distant 
external charges. 

In complete parallel with the KP case one can consider conformal transformations z 
z",e^, etc and map plasma to various geometric configurations. Considering systems with two 
sublattices like ( |30D one arrives at the model of two component plasma in the metal-dielectric 
corner or other bounded regions. This leads again to a boundary Sine-Gordon model at the 
free fermion point corresponding to the temperature F = 2. 



17 



8. BKP hierarchy: a dipole gas on a line between two 
ideal conductors 



In this section we consider real self-similar reductions of the BKP hierarchy corresponding to 
dipole gases on a line between two conductors. 

Let us write out once more the A^-soliton r-function of the BKP hierarchy in the slightly 
different notations 

cr=0,l \^<i<j<N 1=1 / 

^A,. ^ (g^ - aj)ibi - bj){ai - bj){bi - aj) 
(oj + aj)(bi + bj){ai + bj){bi + Oj) ' 

where z,j are the numbers of solitons (they are equivalent in meaning to the coordinate z of 
plasma particles). Choosing 

= exp(— vr/ii/L), 6j = — exp(— 7r(/ii + 'i = l,...,A^, (37) 

we get an expression for the grand partition function of the form ( p!OD for a homogeneous lattice 
gas. The particles of the gas interact via the following VT-potential 

Wd{i - j) = W{h{z - j)) - {h{i - j) -a)-\w {h{i - j) + a) (38) 

with W given by (0). It is seen immediately that this is the potential of two dipole molecules 
consisting of two opposite charges situated at the distance a from each other without the part 
describing interaction of charges inside the dipoles (which is constant). The distance between 
molecules at i-th and j-th site is h{i — j), where h denotes the lattice spacing. The dipoles 
move along the line situated at equal distances 1/2L from two parallel conductors (see the 
Fig. 3). All dipoles are oriented in one direction. Due to the dipole gas interpretation one 
has the effective temperature F = 1 — this is a demonstration of some discrete temperature 
renormalization effect. 

If one takes as in { ^7\) but changes the sign of bi (i.e. shifts + iL), then the signs 

of the second and third terms in (|38[) are changed too. This situation corresponds to a dipole 
gas on the middle line, the dipoles being charged molecules of total charge +2 with the same 
distance between charges a. These molecules are positioned similar to the previous case. 

Another lattice gas model of charged dipoles appears if one replaces a in (|38| ) by i{a+L), < 
a < L. In this case charged dipoles are positioned vertically and symmetrically with respect to 
the middle line y = L/2. The case when a in (^) is replaced by la, < a < L, corresponds to 
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the neutral dipoles gas in the strip between dielectric walls. The dipoles are perpendicular to 
the middle line y = L/2, similar to the previous case, and have an identical orientation. 

The M-periodic self-similar reductions when Oj+M = qai^h+M = Qbi, describe the 

gas consisting of M different types of dipoles. In some particular cases this leads to dipole gases 
with different orientations of neutral or charged molecules in a strip between the conducting 
walls. 

For instance, for the choice 



J2i 



_^-l(2th+a/2) 



e L 



f{{2i+l)h+a/2) 



'2i+l 



-e L 



f{{2i+l)h~a/2) 



neutral dipoles situated on the even and odd sites have opposite directions. (Here one may 
note that the formal substitution CTj (— l)Vj + 1 — (—1)* in the grand partition function 



converts the interaction energy between molecules to the previous form (|38|). However, this 
transformation changes the form of interaction with external fields.) 
The 2M-periodic reduction 



0,2jM-M — 
0'2jM-M+l 



b2jM-M — 
b2jM-M+l 



. _g-f {2j/i+iaM-l) 



n — p-f{2i'i-iai) 

0'2jM-l — e -J I 

Ci2jM 



g-f((2i+l)ft+iai) 



h — p-f{2i/i+iai) 

02jM-l — —e J I 



boi 



•2jM 



.^~i{{2j+l)h^iai) 



g-f ((2i+l)h+iaM-i) 



g-f ((2j+l)/i-iaM-i) 



e L 



f((2j+l)h+ia„) 



b2jM+M-2 - 

02jM+M-l — —e -I ' J^" 



0'2jM+M~2 
0'2jM+M-l 

describes a gas of neutral dipoles lying on the middle line between two ideal dielectrics. Dipoles 
are pointed normally to the boundaries. They can switch their orientations ( "up" and "down" ) 
and internal degrees of freedom characterized by M different dipole moments aj. In general, 
we have to introduce M different chemical potentials describing "internal energy" of the dipole 
molecule. 

Another possible generalization describes mixtures of the +2 charge molecules with both 
parallel and perpendicular orientations of dipoles. Such models can describe polar plasmas 
where molecules can perform discrete rotations by 7r/2. More complicated types of mixtures of 
plasma particles are possible as well. Several physical variables can be calculated here. These 
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Figure 3: BKP: dipole gas on the line between two ideal conductors. The dipoles are neutral, 
oriented identically and lie upon the homogeneous lattice. 

are the number density, polarization and the pressure. One particular model is solved in the 
next section. 

9. Solution of a dipole gas model 

It is known from the theory of solitons that the r-function (^2|) can be represented as a deter- 
minant of some matrix. In general it is necessary to apply the inverse scattering method for an 
auxiliary linear problem to write solutions in such a form. 

Determinant representations are of great help for evaluations of the partition functions, since 
the corresponding matrices appear to have the Toeplitz form in some physically interesting 
cases and, as a result, they can be diagonalized by the discrete Fourier transformation. Here 
we consider only the BKP case, since some models related to the KP hierarchy happen to be 
considered already in the literature, see, e.g. 0. 

As shown in the N = 2p soliton BKP r-function, or the corresponding Ising model par- 
tition function (|33|) admits the following determinant form depending on the soliton parameters 
a{z)^h{z) and the magnetic field (|34[) : 




{a{z) + a{z')){h{z) + Kz')){a{z) - b{z'))ib{z) - a{z')) 
{a{z) - a{z')){b{z) - b{z')){a{z) + b{z')){b{z) + a{z')) 



) 



1/8 



Vdet Gr, 



(39) 



20 



where matrix elements of the matrix Gr are 



b{z) + o(r ) b{z) — a[z') 

^ ^^'^^a(z)-K^') ^ a(z) + a(z') 



and 



n 



\ 1/4 

(a(^) - a{z')){b{z) + b{z')){b{z) - a{z')){a{z) + b{z')) ^ 



^ (aiz) + aiz'))ibiz) - b{z')){biz) + a{z')){a{z) - b{z')) 



(40) 



We consider the situation when the external potential is constant, w = fi, i.e. the dipoles 
interact only between themselves. 

Since the potential (|38|) is of the short range type, the magnetic field in the corresponding 
Ising model ( p4|) is homogeneous in the bulk. We can neglect the field inhomogeneities at the 
edges in the thermodynamic limit oo. /^From (|34|) it follows that 



1 1 °° 

if = -(/i + c), c = -Y.w,{t). 

Consider the dipole gas model described in the Fig. 3. Substituting ( P?] ) into ( pH ) we see 
that the matrix Grjj has the Toeplitz form, Grjj = Grj_j, in the thermodynamic limit ^ oo. 
Therefore it is diagonalized by the discrete Fourier transformation. The final answer for the 
thermodynamic potential per site or the pressure in appropriate units -P(/i) = limp^oo(lnG)/2]9 
can be found from the results of calculations of the free energy per site of the corresponding 
Ising chain 



PU) = - In iQ^Q^bq,q/b; g)oo 

4 {-q,-q,-bq,-q/b;q)^ ' An Jo 



1 

H / du\n\2p{u) 



(41) 



where q = e '^^/^^ b = —e ™/^, and 



p(z/) = cosh(C + /i) + 



2 

oo 



(_e-, -ge— ; g)oo V {b-\qb-q) 
The standard notations for the g-shifted factorials 



(6-le'^ qbe-'^- q)^ ^ (6e'^ gfo-^e"'- q). 



{b,qb i;g), 



(a; q)o = 1, (a; g)„ = (1 - a) (1 - ag) ■ ■ ■ (1 - ag" ) 



(ai, a2, . . . , a„; g)„ = (ai; g)„(a2; g)„ . . . (a^; g)r 
are used in these formulae. 
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Taking the derivative with respect to fi we find the number density of molecules n(/i): 

1 1 ( 1 r du \ 

n(n) = - + - 1 / — — — tanhfC + u), (42) 

2 2V Trio 1 + d(z/) cosh(/x + C) y v p;, v ; 

where 

{qb,q/b;qU{\b\-'/' + \b\'/')e,{u,qy^) 
(-g;g)^2Re^2(i^-(i/2)ln|6|,gV2) " 

Here 62(1^, q^^'^) is the Jacobi theta-function 

00 

e^iu, = 2 J2 cos {2n + l)u 

n=0 

In a similar way one can analyze other dipole gas models, in particular, the general M-periodic 
reduction cases, mentioned in the previous section. 

In conclusion of the discussion of the BKP hierarchy it should be mentioned that statistical 
mechanics interpretation of the KP reductions of the C and D types (CKP and DKP) [1^ is 
not known to the authors. One may try to find a Coulomb gas realization of soliton solutions 
of these and other known integrable equations which were not considered in this paper. 



10. Conclusions 

Constructions considered so far have one essential drawback from the point of view of statistical 
physics: the partition functions derived from the tau functions of classical integrable hierarchies 
correspond only to some fixed temperatures. A possible way of overcoming this obstacle is to 
look for appropriate quantum generalizations of classical hierarchies. 

It is well known (e.g., see [|13|) that the grand partition function of two component plasma 
can be expressed in terms of the following equivalent field theories: Sine-Gordon, Thirring or 
sigma-model. 

Lattice versions of the neutral Coulomb plasma correspond to some discretizations of the 
above models: lattice Sine-Gordon, scalar Hubbard or XXZ model. Since these three models 
are equivalent, we take the scalar Hubbard model as their representative. 

The two component plasma at the temperature ( p2D is mapped onto the free fermion point 
of the one-dimensional scalar Hubbard model. It can also be mapped to free spin 1/2 fermions 
on the lattice 0. The Coulomb plasma at the temperatures different from (p2|) corresponds to 
interacting fermions. 
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/^From our point of view, the relation of Coulomb plasmas to the theory of integrable 
hierarchies described here sheds some new light onto the relation between fermion models and 
different kinds of plasmas. Indeed, it is known [|T^ that the soliton equations can be derived in 
the framework of the free fermion formalism and they are equivalent to Plucker relations on the 
infinite dimensional Grassmanian manifold. We remind briefly basic points of this formalism. 

Let us take fermions on the one-dimensional lattice. Their creation and annihilation oper- 
ators tp*, tpi satisfy the relations: 

The neutral vacuum is deflned as the state where the sites with i < are empty and other sites 
are flUed by fermions (actually, it is a kink state of the XX model or the free fermion point of 
the scalar Hubbard model). 

The action of the fermion operators on the vacuum is 

ipni'va.c) = 0, n < 0, ipnl^^*^) = 0, n >0, 

{vaclipn = 0, n >0, (vac|?/'* =0, n < 0. 
The fermion field operators 

oo oo 

^(A) = y: = E 

fc=— oo fc=— oo 

evolve under the KP flow as follows 

oo 
n=l 

where the generating function of KP Hamiltonians H[t) has the form 

oo oo 

H{t) = Y.H„U, H^= -^k^k+n-, (43) 

n=l k=—oo 

where the colons mean the normal ordering with respect to the vacuum \vac). The r-function 
is deflned as follows 

r{t,g) = (vac|e^W^|vac), g e G, 
where G is a subgroup of the fermion algebra preserving one-fermion states 

G={g\gVg-' = V, gV*g-' =V*} 
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and V = (Bi^zC^pi, V* = (Bi(zzCip*. Choosing the element (7 of G in the definition of r-function 
as 

AT 

g = exp^7iV^(ai)V^*(-6i), 

i=l 

where 7^ are some arbitrary constants related to the zero time phases of solitons, one gets the 
KP iV-soliton r-function (0), (^). 

Thus, the two component plasma at the temperature (p2|) is described simultaneously by 
the KP hierarchy and the free fermions model. Variation of the temperature from this value, 
which in some models (e.g., in the one-dimensional Ising chains picture) is not qualitatively 
distinguished from the other ones, would correspond to a generalization of the free fermion 
formalism for integrable hierarchies to the case of interacting fermions. The KP hamiltonians 
which are of the form of hopping terms of the Hubbard model, suggest that such a generalization 
is possible, provided KP hamiltonians are replaced by conserved quantities of a generalization of 
the spinless Hubbard or XX model. For instance, in the XX limit of the XXZ model derivatives 
of the transfer matrix are in the algebra (|i3| ) 

(lnT(M))(") oc/7_„ + i7„. 

One may conjecture that the derivatives of general XXZ transfer matrix belong to a general- 
ization of the algebra of free fermion Hamiltonians of the KP hierarchy. 

Here we encounter one difficulty. Roughly speaking, the algebra of the KP Hamiltonians 
(|43|) has two times as many Hamiltonians as derivatives of the XXZ transfer matrix. Probably 
there are some extra integrals of motion which are reduced to if„ — in the free fermion 
limit. As an example, we can mention that some extra nonlocal constants of motion have been 



derived for the Hubbard model in |18] 



Another possible way of generalization is due to the approach described, e.g., in [|19|, where 



some generalization of the nonlinear Schrodinger equation is presented. Unfortunately the 
notion of r-function is clearly defined only in the limiting cases of the free fermion points (e.g., 
for the impenetrable Bose gas). It is not obvious that solutions of the corresponding integro- 
differential equations can be expressed in terms of a r-function at general couplings and that 
such a function would make sense in the plasma picture. 

Concluding this article we discuss possible physical significance of the models considered in 
it going beyond the Coulomb gas picture. As was shown in there are nice interpretations of 
the one-dimensional cases from the point of view of Ising magnets. However, various boundary 
conditions arising within the intrinsically two-dimensional Coulomb interaction systems look 
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somewhat artificial in the Ising picture. Still, a number of 2D Ising models with such non-local 
exchange can be formulated which are exactly solvable by the techniques due to Gaudin (see 
the third paper in 0). 

Another possible application concerns the fractional quantum Hall effect (FQHE). Remind 
that the two-dimensional one component Coulomb plasma in the uniform neutralizing back- 
ground of the density l/2m has the following partition function 

Z = J dzi . . . dzN exp {-2mE) , E = -^log\zi - Zj\ + ^Yll^il"^ ■ 

i<j i 

It was shown by Laughlin that the n-point correlation functions in the appropriate state of the 
FQHE at the filling factor 1/m coincide with those of the one component Coulomb plasma. In 
particular, the normalization factor of the corresponding wave function is given by Z. 

Our model is a bit different from the pure one component plasma, since we have non-trivial 
boundary conditions. However, by placing the domain of concentration of the charged particles 
far from the boundary we get the Laughlin plasma. It is also possible that the plasma with 
boundaries has some meaning in the fractional quantum Hall picture as well describing there 
some boundary effects, r-functions of the classical hierarchies correspond to boundary Laughlin 
states at the full filling 1/m = 1. Nontrivial examples m> 1 could be described by r-functions 
of the generalized hierarchies discussed above. 
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